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ABSTRACT 

The  accurate  calculation  of  the  magnetic  fields  steady  or 
rapidly  varying  is  extremely  important  in  designing  pulsed  power 
supplies  for  controlled  thermonuclear  fusion  experiments,  lasers, 
etc.,  where  the  traditional  simplifying  assumptions  become  un¬ 
acceptable  -  especially  when  ferromagnetic  materials  in  high  mag¬ 
netic  fields  are  used. 

A  finite  element  method  for  solution  of  Maxwell's  equations 
for  a  moving  media  in  terms  of  the  magnetic  vector  potential  and 
electrokinetic  scalar  potential  describing  the  penetration  of  the 
magnetic  fields  in  fast  pulsing  power  supplies  of  electromechanical 
type  is  presented.  The  formulation  for  the  steady-state  magnetic 
fields  in  nonlinear  media  results  as  a  particular  case  of  the  me¬ 
thod. 

This  approach  was  used  for  predicting  the  discharge  parameters 
for  the  very  fast  discharging  homopolar  machine  (FDX)  designed  by 
the  Energy  Storage  Group  at  the  University  of  Texas.  FDX  is  in  an 
advanced  state  of  execution. 

This  work  was  supported  by  the  Electric  Power  Research  Insti¬ 
tute  (EPRI)  and  the  Energy  Research  and  Development  Administration 
(ERDA) . 

Introduction 


The  generation  of  powerful  current  pulses  by  pulsed  power  supplies  is 
intimately  related  to  high  magnetic  fields  accompanying  them  and  to  the 
penetration  of  such  currents  and  diffusion  of  such  magnetic  fields  into 
metallic  current  paths  of  the  generator,  transmission  line  and  load. 

An  accurate  prediction  of  the  space  and  time  distribution  of  magnetic 
fields  is  extremely  important  for  the  design  of  pulsed  devices. 

Finite  element  methods  have  been  widely  used  since  1770,  in  solving 
steady  state  problems  of  magnetic  fields,  in  complex,  nonlinear  systems 
[1,2,3]. 

In  such  situations  the  nonlinearity  of  the  ferromagnetic  materials  is 
handled  by  linearizing  the  problem  locally,  solving  the  resulting  linear 
equations  and  introducing  the  nonlinearity  of  the  magnetization  curve  by 
an  iterative  correction  to  the  linearized  equations.  Of  course  this  as¬ 
sumes  that  the  curve  B  =  f (H)  is  single  valued,  neglecting  hysteresis  ef¬ 
fects. 

The  solution  of  the  diffusion  of  magnetic  fields  by  a  finite  element 
method  -  Galerkin  projective  technique  was  attempted  in  several  papers 
[4,5].  Foggia  et  al  [4]  assumed  a  steady  state  regime  and  a  prescribed 
sinusoidal  travelling  field  impressed  to  the  physical  model.  Miya  [5]  has 
a  similar  treatment  prescribing  the  variation  of  plasma  current  in  time. 
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Actually  It  is  not  the  solution  for  the  forced  response  of  the  system  but 
the  solution  for  the  "natural  response"  or  rather  the  "complete  response" 
that  is  required  for  pulsed  power  systems.  The  finite  element  method  pre 
sented  herein  solves  the  transient  diffusion  problem  in  terms  of  a  magne¬ 
tic  vector  potential  A  through  an  intermediary  electrokinetic  scalar  po¬ 
tential  ¥  which  transforms  the  forcing  function  concentrated  in  the  gene¬ 
rator  into  a  distributed  forcing  function  which  acts  locally  in  each  poin 
of  the  system:  generator,  transmission  line,  and  load. 

The  physical  problem 

The  physical  problem  treated  herein  is  the  discharging  of  a  very  fas 
homopolar  generator  which  is  a  variant  of  the  FDX  experimental  machine 
presented  in  detail  in  a  companion  paper  at  this  conference.  The  scope  o 
FDX  is  to  explore  the  fundamental  limitations  of  fast  discharging  homopo¬ 
lar  power  supplies [6]. 

The  machine  consists  of  two  counterrotating  aluminum  rotors,  storing 
at  rated  speed  a  kinetic  energy  of  0.36  Mi. 


Figure  1:  Machine  configuration 


The  excitation  is  provided  by  a  four  turns  field  coil  pulsed  by  discharg¬ 
ing  a  5  Mj  slow  discharge  homopolar  generator.  The  duration  of  the  exci¬ 
tation  pulse  is  long  compared  with  the  time  of  discharge  of  the  FDX  ma¬ 
chine  (the  pulse  has  almost  a  flat  top  for  120  milliseconds).  For  this 
reason  the  flux  density  of  the  excitation  field  is  assumed  to  be  constanl 


IIIE2-3 


in  time.  The  flux  linking  the  active  area  of  a  rotor  is  0.218  Webers  at  an 
average  flux  density  of  4  Tesla  which  corresponds  to  a  no-load  voltage  of 
104  V  per  rotor,  208  V  per  machine  calculated  as  line  integral  of  V  x  BQ. 


Figure  2:  Excitation  flux  lines  in  rotor 


Figure  2  shows  a  map  of  excitation  flux  linking  the  rotor,  at  the  mo¬ 
ment  when  the  peak  value  is  reached.  At  this  time  the  brushes  are  applied 
to  the  rotor  and  a  moment  later  the  making  switch  is  closed.  The  dis¬ 
charge  current  begins  to  flow,  at  the  beginning,  in  a  very  thin  layer  or 
rotor  and  compensation  plates,  moments  later  penetrating  into  the  interior 
of  conductors. 

The  J  x  Bq  forces  stop  the  rotor,  accelerate  it  in  the  opposite  di¬ 
rection,  stop  it  again  until  the  stored  energy  is  dissipated  in  space  and 
time  of  the  electromagnetic  field  in  the  process  described  above  is  the 
subject  of  this  work. 

The  mathematical  electromagnetic  problem 

The  time  scale  of  the  electromagnetic  phenomena  concerning  this  par¬ 
ticular  problem  the  assumption  of  a  quasistationary  state  holds  and  Max¬ 
well's  equations  are  used  accordingly. 

The  energy  stored  at  t  =  0  in  different  parts  of  the  system  expressed 
through  the  electromagnetic  state  variable  represents  the  input  for  Max¬ 
well's  equations,  and  the  "natural  response"  of  this  system  is  described 
by  the  solution  sought  through  the  finite  element  approach. 

The  aforementioned  treatments  of  the  diffusion  of  electromagnetic 
field  [4,5]  were  seeking  the  "forced  response"  of  the  system.  Boundary 
conditions  are  given  at  all  times  in  such  cases,  and  through  them,  for 
given  material  properties  of  the  system  -  the  electromagnetic  field  quan¬ 
tities  are  uniquely  defined. 

Maxwell's  equations,  for  a  quasistationary  regime,  can  be  written  as 


(a)  V  x  h  =  J 
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3B 

(b) 

V  x  E  =  - 

"  3t 

(] 

(c) 

VB  =  0 

Ohm' s  law  in  a  moving  media  is  given  by 

J  =  cr(E  +  V  x  b)  (: 

and,  in  view  of  the  assumption  of  homogeneous,  isotropic  and  hysteresis 
free  materials 


B  =  pH  ( 

In  order  to  distribute  the  effects_of  _the  stored  energy  in  the  rotor  cha 
acterized  by  the  line  integral  of  V  x  B  as  a  variable  of  state,  we  intro¬ 
duce  the  electrokinetic  scalar  potential  ¥  as  an  intermediate  variable. 
In  the  moment  of  closing  the  switch,  the  electrokinetic  state,  character¬ 
ized  by  ¥  is  transmitted  (with  the  speed  of  light)  to  each  point  of  the 
closed  electrical  circuit. 

This  choice  forces  us  to  select  as  primary  independent  variable  the 
magnetic  vector  potential  A,  leading  to 

9A 

!  '  -V  '  at  < 

B  -  V  x  A  ( 

To  these  equations  we  adjoin  the  conditions  that  div  J  =  0,  which  expres 
es  the  conservation  of  charge  condition  and  div  A  =  0. 

Substitution  and  combining  the  above  equations  yields  the  set 

OA+Vx-jjV  x  A  -  a(V  x  V  xa  -  VT)  =  0 

7'OV?  -  V*a(V  xVxA-A)=0  ( 

subject  to 


V* A  =  0. 

Because  both  the  rotor  motion  V  and  the  potential  of  the  applied  ma 
netic  field  Aq  are  axisymmetric  and  have  only  9-components,  the  problem 
separates  into  two  weakly  coupled  equations,  (7)  which  governs  the  distr 
bution  of  the  steadily  applied  field  and  (6)  in  which  the  vector  A  has  o 
ly  r  and  z  components. 

Vx-Vxa  =J  ( 

~  y  ~  ~o  ~o 

The  finite  element  solution  of  (7),  although  necessary  to  our  probl 
is  a  standard  one  and  need  not  be  discussed  further  here. 

In  addition  to  the  equations  (6),  we  have  as  side  conditions,  the  v 
nishing  of  the  potentials  at  large  distances  from  the  conductors  and  the 
usual  continuity  or  jump  conditions  on  current  and  field.  Because  the 
current  distributions  do  not  extend  to  infinity: 
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and 


A  =  A  =  Y  =  0  as  r  00 
r  z 


(8) 


n  x 


T7<V  x  A) 

U  - 


-  0 


n 


a(v  x  v  X  A  -  VY  -  A) 
~o 


=  0 


(9) 


In  (9)  the  brackets  denote  the  jump  in  the  enclosed  quantities  across  an 
interface  between  different  materials  or  between  moving  and  stationary  ma¬ 
terials.  In  the  finite  element  model  of  the  region  shown  in  Figure  1,  the 
symmetry  about  r  =  0  and  the  antisymmetry  about  z  =  0  give  rise  to  the 
following  additional  boundary  conditions 


A  =  0 
r 


on  r  =  0 


A  =  Y  =  0  on  z  =  0, 
r 


Finite  Element  Formulation 


(10) 


We  seek  solutions  to  (6)  and  (9)  subject  to  the  conditions  (8)  and 
(10)  in  the  weak  or  Galerkin  sense.  To  this  end  we  introduce  a  Lagrange 
multiplier  X  corresponding  to  the  constraint  (6c)  and  the  test  functions 
u,  r),  and  T  corresponding  to  A,  Y,  and  X.  The  weak  statement  of  the  prob¬ 
lem  is 

Z{//[cru*A  +  —  (V  x  a)  •  (V  x  u)  +  au*VY  +  XV*u  +  tV*A 

i  fl.  '  ~  -  U  ~  ~  ~~ 

1  -a u*(V  X  v  X  A  )]da}  =  0 

~  ~  ~  ~o 

and  (11) 


£{//[ctVtvVY  +  crA*Vri  -  a(V  x  v  x  a  )*Vr|]da}  =  0 
i  fl.  ~  -  ~  ~~~o 

x 

for  arbitrary  u,  n>  and  T  vanishing  at  °°. 

It  is  significant  that  the  jump  terms  have  equilibrated  the  boundary 
integrals  from  fhe  integration  by  parts,  i.e.  the  jump  conditions  are 
natural  conditions  that  follow  from  (11) . 

The  discrete  form  of  (11)  is  obtained  by  introducing  the  finite  ele¬ 
ment  approximation  of  the  fields  A,  u,  Y,  n,  X  and  T.  The  elements  are 
quadrilaterals  and  the  shape  functions  are  quadratic  for  A,  u,  Y,  and  n  and 
linear  for  X  and  T.  Thus  each  comer  node  of  the  mesh  hai  a~total  of  four 
unknowns,  i.e.  A  ,  A  ,  Y  and  X;  while  the  midside  nodes  have  only  A  ,  A  , 
and  Y.  If  the  set  of  nodal  point  values  is  denoted  q  then  the  discrete2 
form  of  (11)  can  be  written  as 

C  q  +  K  q  =  f  (12) 


In  (12)  C  and  K  are,  for  an  element,  28  x  28  matrices  whose  elements  are 
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the  integrals  of  products  of  shape  functions.  It  is  apparent  from  (11) 
that  neither  £  nor  g  is  symmetric.  The  forcing  function  f  in  (12)  con¬ 
tains  terms  proportional  to  V  x  V  x  A  .  The  motion  of  the  rotor  ia  treat¬ 
ed  as  given  in  this  analysis'"and~the~value  of  the  applied  field  A  is 
known  from  an  initial  finite  element  solution  for  the  single  compSttent 


Results 


The  first  step  in  the  analysis  is  to  calculate  steady  excitation 
field  B  in  which  the  discharge  takes  place.  The  forcing  function  respons 
ible  for  discharge  (energy  conversion  arid  transfer) ,  [V  *  B  ]  has  radial 
and  axial  components  corresponding  to  the  so-called  "dlsk"~and  "drum"  con¬ 
figurations  of  homopolar  machines. 

The  z  component  is  very  small  when  compared  to  the  radial  one  but  in 
the  area  covered  by  brushes  it  can  cause  important  circulating  currents 
along  the  rotor-brush  interface  (Figures  3,  4). 


Lgure  3: Magnetic  field  BQ  t  =  5  ysec. 

u 

Figures  3  and  4  show  the  contours  of  Bq  at  an  early  and  late  time  during 
discharge.  The  current  density  vector  J  with  no  component  in  the  0  direc¬ 
tion  produces  a  flux  density  with  only  one  component.  J  is  tangent  to  the 
contours  and  its  magnitude  is  inversely  proportional  to  the  distance  be¬ 
tween  them. 

At  5  ysec.  all  the  current,  hence  all  the  field,  is  confined  to  the 
first  layer  of  elements  adjacent  to  the  airgap.  Later  (Figure  4)  the  cur¬ 
rent  and  field  penetrate  gradually  -  tending  toward  a  uniform  distributioi 
At  the  final  time  the  rotor  has  stopped  and  begins  t^  reverse.  The 
J  x  B(N/m^),  electromagnetic  force  densities,  and  pJ  (W/m^),  specific  ohm: 
loss,  are  easily  available  from  the  field  distributions.  From  them  desigi 
parameters  as  stresses  and  heating  as  functions  of  time  can  be  calculated. 
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Figure  4:  Magnetic  field  Bfl  t  =  250  ysec. 
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